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ABSTRACT 

We show that there exists a maximal possible accretion rate, beyond which 
global slim disk solutions cannot be constructed because in the vertical direction 
the gravitational force would be unable to balance the pressure force to gather 
the accreted matter. The principle for this restriction is the same as that for the 
Eddington luminosity and the corresponding critical accretion rate, which were 
derived for spherical accretion by considering the same force balance in the radial 
direction. If the assumption of hydrostatic equilibrium is waived and vertical 
motion is included, this restriction may become even more serious as the value of 
the maximal possible accretion rate becomes smaller. Previous understanding in 
the literature that global slim disk solutions could stand for any large accretion 
rates is due to the overestimation of the vertical gravitational force by using an 
approximate potential. For accretion flows with large accretion rates at large 
radii, outflows seem unavoidable in order for the accretion flow to reduce the 
accretion rate and follow a global solution till the central black hole. 

Subject headings: accretion, accretion disks - black hole physics - hydrodynamics 

1. INTRODUCTION 

Although great progress has been made in recent years in increasingly sophisticated 
numerical accretion disk simulations, simple analytic disk models still are the only accessible 
way of making direct link between the theory and observations, as only these models can 
be used to estimate, e.g., the spectra of accretion-powered astrophysical systems. This fact 
justifies the continuous effort to improve the understanding of accretion processes using a 
simple viscosity parameterization. 

The slim disk model is one of such simple analytic models for black hole accretion 
(Abramowicz et al. 1988; Kato et al. 1998). This model was developed upon the standard 
Shakura-Sunyaev disk (SSD) model (Shakura & Sunyaev 1973) by considering two processes 
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in black hole accretion flows, namely the transonic motion and the advective heat transport, 
which were neglected in the SSD model. In the limit of low mass accretion rates, i.e., M is 
substantially lower than its critical value MEdd corresponding to the Eddington luminosity 
Lgdd) say, M < O.lMEdd) the advective heat transport is unimportant, and the structure of 
slim disks is similar to that of SSDs, with a difference that slim disk flows are transonic in 
their inner regions (e.g., Chen & Taam 1993). In this sense, the theoretical basis of black 
hole disk models for low and even moderate accretion rates is well established. But the slim 
disk model is supposed to have the advantage over the SSD model that it can be extended 
to the case of high accretion rates, i.e., with M approaching or surpassing Medd- In this 
case, the very basic assumption of the SSD model that is vahd only for low accretion rates, 
i.e., the geometrical thinness, H r, where H is the half-thickncss of the disk and r is 
the cylindrical radius, would break down (Frank et al. 2002, p. 98). It was then suggested 
in the slim disk model that the disk becomes geometrically slim, i.e., with H < r; and 
accordingly, the process of advective heat transport becomes important or even dominant 
over the radiative cooling because, as given by Abramowicz et al. (1986), there is a relation 
/adv « {H/rf, where /adv = QZiv/Qtis is the advective factor, with Q^^ and being the 
advective cooling and viscous heating rates per unit area of the disk, respectively. 

However, the self-consistency of the slim disk model and its applicability to the high 
accretion rate case are not so obvious. In particular, uncertainties seem to concentrate on 
the treatment of the vertical structure of the disk. First, vertical hydrostatic equilibrium 
is a reasonable assumption for SSDs because, for geometrically thin disks, vertical motion 
of the disk matter must be negligible compared with radial motion; but it is questionable 
to adopt this assumption for slim disks that are not thin. Second, as noticed recently by 
Gu & Lu (2007, hereafter GL07), even the vertical hydrostatic equilibrium is assumed, in 
the slim disk model there is a serious inconsistency that the vertical gravitational force was 
greatly magnified by using the approximate form of potential due to Hoshi (1977), which is 
valid only for thin disks again; and accordingly, the previous understanding that slim disks 
could exist for any large accretion rates seems doubtful. Third, the above mentioned relation 
/adv {H/tY was derived with the Hoshi approximation of potential, so its applicability 
to not thin disks with large accretion rates is not justified; or in other words, it is not clear 
whether an accretion disk can ensure advection dominance (i.e., /adv > 0.5) while remaining 
to be geometrically slim. According to the analyses of Narayan & Yi (1995) and Gu et al. 
(2008), in order for advection to be dominant, the disk needs to be geometrically thick, i.e., 
with H > r, rather than slim. This is probably the reason why H > r was obtained in many 
numerical calculations of slim disks and other accretion disks (see references in GL07). 

What about another popular black hole accretion disk model, the advection-dominated 
accretion flow (ADAF) model (Narayan & Yi 1994; Abramowicz et al. 1995)? ADAFs are 
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also supposed to be geometrically slim and advection-dominated, and their vertical struc- 
ture was treated in a way similar to that for slim disks. The above mentioned uncertain- 
ties should also apply to the ADAF model, because they are based on purely hydrostatic 
considerations and are related only to the geometrical thickness of the disk. Many two- 
and three-dimensional numerical simulations of viscous radiatively inefficient accretion flows 
revealed the existence of convection-dominated accretion flows rather than ADAFs (see ref- 
erences in GL07). This fact is probably an indication that the ADAF model might have 
hidden inconsistencies, and one of which might be related to the treatment of vertical struc- 
ture. What is different from slim disks is that ADAFs are expected to correspond to very 
low accretion rates and are known to have a maximal possible accretion rate at each radius 
(e.g., Abramowicz et al. 1995), so the problem addressed in GL07 regarding the allowed 
accretion rate may have no impact on ADAFs. In addition, ADAFs are optically thin and 
ion pressure-supported, the radiation processes in them are more complicated, making the 
vertical structure more difficult to deal with, than for slim disks that are optically thick and 
radiation pressure-supported. 

Our present work is devoted to discuss the slim disk model and is a straightforward 
continuation of GL07. All the results of GL07 were based on a local analysis, i.e., only for 
a certain radius. Although a similar local analysis was often used in the literature (e.g., 
Abramowicz et al. 1995; Chen et al. 1995; Kato et al. 1998), one should be cautious of the 
fact that conclusions made in the local sense do not necessarily hold in the global sense. For 
example, it has been shown that, even though a disk is locally unstable at a given radius 
according to a local stabihty analysis, the disk can be globally stable in global numerical 
simulations (Janiuk et al. 2002; Gicrlihski & Done 2004). Therefore, it is worthwhile to 
check and extend the results of GL07 by investigating global solutions of original differential 
equations for black hole accretion flows, similar to what was done by, e.g., Chen & Wang 
(2004), Watarai et al. (2005), Artemova et al. (2006), and Watarai (2006), but with a 
revised vertical gravitational force. 

2. SOLUTIONS WITH VERTICAL HYDROSTATIC EQUILIBRIUM 

2.1. Equations 

The basic equations to be solved for shm disks can be written in cyhndrical coordinates 
as (cf. Kato et al. 1998, p. 236; GL07): 



M — —2'KrTiVr — constant. 



(1) 



-4- 



^.^ + -^ + (^l|-^^> = 0, (2) 

dr po dr 

M{nr^ - j) = 27iar^U, (3) 

Q\is — Qadv + ^rad' i^) 

Po = \- o > (5) 

where E = 2 jj^ pd^; is the surface density, p is the density, p is the pressure, Vr is the radial 
velocity, is the angular velocity and Qk is its Keplerian value, j is an integration constant 
representing the specific angular momentum accreted by the black hole, a is the Shakura- 
Sunyaev viscosity parameter, H = 2 pdz is the vertically integrated pressure, is the 
radiative cooling rate per unit area, T is the temperature, p, is the mean molecular weight 
and is taken to be 0.62, and the subscript "0" represents quantities on the equatorial plane. 

There are some differences between equations (1-5) and the basic equations in GL07 
(their eqs. [6 - 10]). In GL07 all the five equations were written in the vertically integrated 
form; here the continuity equations (1), angular momentum equation (3), and energy equa- 
tion (4) are in the same form, because it is obviously convenient to write equation (1) in 
this form and equations (3) and (4) also contain M; but the radial momentum equation (2) 
and state equation (5) are given for equatorial quantities. The state equation of GL07 (their 
eq. [10]) is a trivial vertical integration of our equation (5) and makes no difference. But in 
GL07 (and also in Kato et al. 1998), in obtaining the vertically integrated radial momentum 
equation (eq. [7] of GL07 or eq. [8.25] of Kato et al. 1998), some simphfications were made. 
First, both Vr and were regarded to be independent of the coordinate z, and this is also 
adopted in our equations here. Second, a step — W taken, but in the strict 

sense this equality holds only when H does not vary with r, in studies of global solutions 
here we should better not to take it. This is the reason why we prefer to write equation (2) 
for equatorial quantities. The differential equation (4) was reduced to be an algebraic one 
in GL07; and here instead, it keeps its original form with the explicit expressions for Q^^^, 
<5adv' and Q;^^ given as 

_ _Mn{nr^ - j) dlnn 
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n- 3M n / c?lnpo A d\npo . 

^^"^ 27rr2E \ dhir 3 dlnr ' ' ^ ' 

(Kato et al. 1998), where k = k^^ + ks = 0.34 + 6.4 x lO^^pT "^^^cm^ g-i, p(= E/2i/) and 
T(= Tdz/H) are the vertically averaged density and temperature, respectively. 

The key difference between the existing slim disk model and our work here is in the 
treatment of the vertical hydrostatic equilibrium equation, 

dp dip 

where i/j is the potential. In the shm disk model (e.g., Kato et al. 1998, p. 241), the potential 
of Paczyhski & Wiita (1980), 

. , GM , , 

Vr"^ ^ — Tg 

where Vg = IGMjc? is the gravitational radius, was approximated in the form of Hoshi 
(1977), i.e., 

^(r,^)^^(r,0) + ^. (11) 

Using equation (11) and assuming a polytropic relation in the vertical direction, p = Kp^^^l^ ^ 
where K and N are constants, the vertical integration of equation (9) gave (for N — 

\ 1/3 / N 1/4 rr 2 

L) ^(L) = (12) 



and E and 11 were 



8^ = ^llH\ (13) 
Po 



E = 2 X ^^poH, (14) 
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n = 2 X ^PoH. (15) 

The simple relation Cs/^kH — constant was obtained, with the sound speed Cg being defined 
either as = Po/po or as = II/E. 

However, as the main point made in GL07, the approximation of equation (11) is invalid 
for slim disks. By using the explicit Paczyhski & Wiita potential, equation (10), to integrate 
equation (9), one obtains instead of equations (12) and (13): 

/ \ 1/3 / \ 1/4 m 1 . . 1 . 

9=[-) = - = 1 1 (16) 



4^ = GM L_ ) . (17) 



Accordingly, E and 11 become 



H 

3 



2po / g'dz, (18) 



H 

4 



n = 2po / g^'dz, (19) 



which are much more complicated than equations (14) and (15). It is seen that the relation 
Cs/^kH = constant does not hold. 

To summarize, in the existing slim disk model, the set of eight equations, namely equa- 
tions (1 - 5) and (13 - 15), could be solved for the eight unknown quantities po, Po, S, IT, Tq, 
Vri f^, and H as functions of r, with given constant parameters M, M, a, and j; while in 
our work here, equations (17 - 19), instead of equations (13 - 15), still along with equations 
(1-5) form a new set to be solved for the same eight unknowns. There are two differential 
equations, i.e., equations (2) and (4); and other equations in each set are algebraic, including 
equations (18) and (19) once the vertical integration is analytically made. 



2.2. Solutions 



We numerically solve the set of equations (1 - 5) and (17 - 19), i.e., with the correctly 
calculated vertical gravitational force using the explicit Paczyhski & Wiita potential. For 
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comparisons, we also solve the set of equations (1 - 5) and (13 - 15), as done already in 
the existing slim disk model, i.e., with the magnified vertical gravitational force using the 
Hoshi approximation of the Paczyhski & Wiita potential. After some algebra, the equations 
in each set can be combined into only two differential equations. A physically acceptable 
global solution should be able to extend from a large radius to the vicinity of the central 
black hole, passing the sonic point regularly. We start the integration at r = lO^Vg, where the 
outer boundary conditions are set to be corresponding to an SSD, i.e., being geometrically 
thin and Keplerian rotating. If a solution is found to be transonic, we extend it to r = 2rg 
(in the following figures only a smaller radial range is shown in order to see more clearly 
the solution behavior in the inner region). We fix M = IOMq and a — 0.1, and vary the 
values of M for different solutions. For a transonic solution, the other constant parameter j 
is an eigenvalue of the problem and has to be adjusted correctly and accurately (e.g., Chen 
& Taam 1993; Chen & Wang 2004). 

Rather than showing the radial variation of every physical quantity, our main purpose 
here is to see how the correction of the vertical gravitational force changes the previous 
understanding of global slim disk solutions, that is, there was no upper limit of M for slim 
disks, any large value of M could correspond to a global solution (e.g., Kato et al. 1998; Chen 
& Wang 2004; Watarai et al. 2005; Watarai 2006). Figures 1 and 2 are sufficient for this 
purpose, which are for the relative thickness H/r and the advective factor /adv, respectively. 
In the figures, solid lines are solutions of equations (1-5) and (17 - 19), and dashed fines are 
solutions of equations (1 - 5) and (13 - 15). Sofid line a and dashed fine a' are for the same 
accretion rate m—1, where m = M /M-^^^, with MEdd = Mt^GM / ck^s being the Eddington 
accretion rate; and similarly, solid line h and dashed line h\ solid line c and dashed line c', 
solid line d and dashed line d\ and solid line e and dashed line e' arc for rh = 10, 31.5, 31.7. 
and 100, respectively. It is seen that for a moderate accretion rate m = 1 (lines a and 
a'), global solutions obtained from the two sets of equations are almost the same, i.e., both 
the solutions have H/r < 1 (Fig. 1), and arc radiation-dominated in the outer region and 
with important advection in the inner region (Fig. 2). This proves that using the Hoshi 
approximation of potential for moderate accretion rates is acceptable. However, deviations 
appear and become serious as rh increases. For rh — 10, our new solution (lines b in the two 
figures) has H and /adv significantly larger than that in the existing slim disk solution (lines 
6'). Such a situation continues till a critical value rh — 31.5, for which our new solution 
can still be constructed (lines c). For a slightly larger m = 31.7 and a still larger m = 100, 
global solutions can be no longer obtained from our new set of equations (1 - 5) and (17 - 
19). It is seen that H tends to infinity (lines d and e in Fig. 1) and /adv tends to exceed 1 
(its maximal possible value, lines d and e in Fig. 2), so the inward integration cannot go on. 
These solutions (fines d and e) are not global solutions at all, since they cannot extend to a 
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sonic point. On the other hand, no matter how large M is, global solutions can always be 
found from the set of existing slim disk equations (1 - 5) and (13 - 15), as drawn by lines b\ 
c', d\ and e' in the two figures (lines c' and d' coincide with each other). 

Unfortunately, the existing global slim disk solutions, though being formally constructed 
for any high accretion rates, have hidden inconsistencies. As seen clearly in Figure 1 of GL07, 
the vertical gravitational force, dip/dz in equation (9), was greatly overestimated, and ac- 
cordingly, the geometrical thickness H was greatly underestimated, by using the approxi- 
mation of equation (11). This is the reason why in the slim disk model, the gravitational 
force seemed to be always able to balance the pressure force and ensure vertical hydrostatic 
equilibrium. Even so, slim disk solutions may still have H/r > 1 (lines c', d\ and e' in Fig. 
1), the assumption of slimness is violated. 

As the main result of our work, it is found that, when the vertical gravitational force 
is correctly calculated from equation (10), there exists a maximal possible accretion rate, 
Mnax ~ 31.5MEdd, bcyoud which there are no global solutions at all (lines d and e in Figs. 
1 and 2). The physical reason for this is the following. The amount of accreted matter that 
can be gathered by the black hole's gravitational force must be hmited. If M > Mmax, the 
pressure force of the matter and radiation in the vertical direction would be too large to be 
balanced by the gravitational force, the disk would be huffed by the pressure force to tend 
to an infinite thickness (lines d and e in Fig. 1), and the accretion processes would never be 
maintained. This reason is based on the same principle as that for the Eddington luminosity 
Ledd and the corresponding critical accretion rate Mgdd to be defined. The only difference 
is that i^Edd and M^dd were derived for spherical accretion, i.e., by considering the balance 
between the gravitational force and the pressure force in the radial direction; and here the 
maximal possible accretion rate is found by considering the balance of the same two forces 
in the vertical direction of accretion disks. 

In GL07, by a local analysis, i.e., considering the vertical balance of gravitational and 
pressure forces at a certain radius, a similar maximal possible accretion rate, Minax('"), was 
found for each radius. Their result is confirmed here. But a global solution is with a 
constant accretion rate, so the allowed accretion rate has a unique value, rather than being 
radius-dependent . 

As a check of our global solutions, the total optical depth r = is calculated and 

is shown in Figure 3, where lines a, b, c, a\ b\ and c' are the correspondents of lines a, 6, 
c, a', b\ and c' in Figures 1 and 2, respectively. It is seen that, similar to the existing slim 
disk solutions (lines a', h\ and c'), our global solutions are also optically thick everywhere 
(lines a, 6, and c). But r in our solutions is somewhat smaller than that in the existing 
slim disk solutions, especially in the inner regions. This is because in our solutions, the 
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correctly calculated vertical gravitational force is smaller, then H is larger, Vr is larger, and 
E is smaller for the same M; and k is almost unchanged, especially in the inner regions 
where the electron scattering opacity is dominant and is a constant. 



3. DISCUSSION 

3.1. Waiving hydrostatic equilibrium 

Though global solutions with vertical gravitational force correctly calculated can be 
obtained for M < M^ax as shown in Figures 1, 2, and 3, there seems to be also some 
inconsistency hidden in these solutions. For high accretion rates, the disk's relative thickness 
H/r becomes substantially larger than 1 (lines b and c in Fig. 1); and for the critical value 
M — Sl.SMEdd; H/r reaches up to ~ 20 in the middle region of the solution. As argued 
by Abramowicz et al. (1997), if it is assumed that there is no velocity component in the 
direction orthogonal to the surface of the disk (i.e., no outflows leaving the disk), then there 
is a relation vn/vr = dH/dr, where vh is the vertical velocity on the surface, and vertical 
motion should not be neglected if H does not very slowly vary with r. In particular, at 
the maximum of H/r (see lines b and c in Fig. 1), it follows from d{H/r)/dr = that 
dH/dr = H/r. This means that in the middle region of the solution the vertical velocity 
could greatly exceed the radial velocity, and the assumption of hydrostatic equilibrium would 
not be vahd. In this case, instead of hydrostatic equilibrium equation (9), one should use 
the more general form of vertical momentum equation. 



1 dp dih dvz dvz , . 

+ + ^'■^ + = 20 

paz oz or oz 



(Abramowicz et al. 1997). 



Solving this partial differential equation is beyond the capacity of the slim disk model 
that is one-dimensional. In an illustrative sense, here we wish to try to consider the non- 
negligible vertical velocity Vz with a very simple treatment, which is similar to what was 
done in Abramowicz et al. (1997). We assume 



where 



z z 

Vz{r, z) = —Vh = -u, (21) 



, , dhiH 

„(r)=..— , (22) 
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then equation (20) is reduced to be 



where 

If du 



Y(r) = ^\^rvr^-vru + u'j, (24) 

which looks similar to equation (9) of Abramowicz et al. (1997). 

Instead of equations (16) and (17), the vertical integration of equation (23) gives (of 
course, still with the explicit potential eq. [10]) 



A^ = GM[ , + -YH\ (26) 

Equations (18) and (19) are formally unchanged, but the quantity g in these two equations 
is given by equation (25) now. 

The nine equations, i.e., equations (1 - 5), (18), (19), (22), and (26) can be solved for 
nine unknowns, po, Po, S, 11, Tq, Vr, ^, H and u (eq. [1] keeps unchanged even there is 
a non-zero v^, because of the no outflow assumption; and eqs. [2 - 5] still hold as well). 
The situation of solutions is seen in Figure 4 that shows H/r as a function of r. Global 
solutions exist till a new maximal possible accretion rate Mmax ~ 8.5MEdd (lines a and b 
that are for m = 1 and 8.5, respectively), above which there are no global solutions (lines 
c and d that are for m = 9 and 30, respectively). But this time, for M > Mmax, it is not 
the case of Figure 1 that the inward integration stops at a radius where the thickness H 
tends to infinity; rather, it is the case that the inward integration stops at a radius where the 
quantity g of equation (25) tends to become unphysically negative at some height z between 
the equatorial plane z = (where g = 1) and the disk surface z = H (where g = 0), even 
though H is finite at that radius. It is not surprising that the value of Mmax here is even 
smaller than that in the above section where vertical hydrostatic equilibrium is assumed. 
The reason for this is the following. As seen from Figure 5, in global solutions a and 6, vh 
is always negative, i.e., the vertical motion is always inward towards the equatorial plane; 
and the absolute value of vh increases with decreasing r. This means clearly that there 
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ls a vertical acceleration towards the equatorial plane. That is, in the vertical direction, 
the gravitational force has to overcome the pressure force to accelerate the accreted matter, 
rather than only balancing the pressure force to ensure hydrostatic equilibrium. The upper 
limit of the amount of the accreted matter, to which the gravitational force is able to do this 
harder job, must be smaller than in the hydrostatic equilibrium case. For M > Mmax, the 
gravitational force would be unable to do the job. To keep equation (23) formally holding, 
g (and p and T) of equation (25) would have to become mathematically negative, which are 
physically unacceptable. 

Though the treatment of vertical motion made here is crude, it is sayablc that inclusion 
of this motion should not change the main conclusion reached in the above section, i.e., there 
is a maximal possible accretion rate Mmax for global slim disk solutions to exist; and the 
value of Mmax becomes even smaller than in the case that this motion was omitted. 



3.2. Outflows 

If the mass accretion rate given at large radii exceeds its maximal possible value, M > 
Mmax, it seems that the only possible way for a global slim disk solution to be realized is 
that the accretion flow loses its matter at large radii in the form of outflows, such that at 
smaller radii M is reduced to be below Mmax- 

Outflows have been observed in many high energy astrophysical systems that are be- 
lieved to be powered by black hole accretion, but the mechanism of outflow formation remains 
unclear from the theoretical point of view. At the same time when the ADAF model was 
proposed, it was suggested that ADAFs are likely to be able to produce outflows because 
they have a positive BernouUi constant in their self-similar solutions (Narayan & Yi 1994). 
Later, Abramowicz et al. (2000) showed that in global solutions of ADAFs the Bernoulli 
function, rather than the Bernoulli constant, can have either positive or negative values, and 
that even a positive Bernoulli function is only a necessary, not a sufficient, condition for 
outflow formation. 

Historically, the Bernoulli constant or function was defined as the sum of the specific 
enthalpy, kinetic energy, and gravitational potential energy. For ADAFs Narayan & Yi (1994) 
and Abramowicz et al. (2000) discussed, the Bernoulli function has a simple expression 
since the enthalpy is totally due to the disk gas; it has also a clear physical meaning that 
its positivity may imply a possibility of outflows. For slim disks, however, the contribution 
from radiation to the enthalpy becomes important or even dominant over that from gas. 
To our knowledge, in this case it is unclear whether and how can the Bernoulli function be 
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defined, or whether this function would have the same meaning as for ADAFs if its historical 
definition is copied. Therefore, we use instead the specific total energy of the matter in slim 
disks to consider the possibility of outflows, whose equatorial-plane value is 

\2 firup po J 2 ^ r -rg 

(cf. eq. [11.33] of Kato et al. 1998). 

Figure 6 shows E corresponding to the solutions with the correct gravitational force in 
Figures 1 and 2. Lines a, b, c, d, and e in Figure 6 are for the same physical parameters as 
lines a, b, c, d, and e in Figures 1 and 2, respectively. For accretion rates that allow global 
solutions to exist, M = 1,10 and 31.5MEdd (lines a, b. and c), E is negative everywhere 
in the disk, and outflows arc unlikely to originate. However, for accretion rates exceeding 
Mmax('"), M = 31.7 and lOOMEdd (lines d and e), it is seen that as r decreases, E tends to 
become positive flrst, and then the solution stops extending. Further, it is noticeable that, 
when a flow changes from a state represented by line e to a state represented by line d, its E 
increases, and at the same time its M decreases. The increase of E favors outflow formation, 
and the decrease of M is just the result of outflows. Such a process of the change of the 
flow's state continues until M is sufficiently reduced and the fiow reaches a state represented 
by line c, then E becomes negative and outfiows cease. 

To summarize, in this subsection we make two arguments for outflows produced in an 

accretion flow with large M at large radii. First, outflows are necessary because the accretion 
flow must lose its matter in this way in order to follow a global solution till the central black 
hole. Second, outflows are possible because the accretion flow can have a positive specific 
total energy; what is additionally needed for outflows to be realized is that the matter in the 
accretion flow comes under an outward perturbation in the vertical direction. 
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Fig. 1. — Disk's relative thickness H/r as a function of r. Solid lines represent our results 
with the vertical gravitational force revised, and dashed lines the results in the existing slim 
disk model. Lines a and a\ b and b\ c and c', d and d\ and e and e' are for accretion rates 
M/MEdd = 1, 10, 31.5, 31.7, and 100, respectively. There exists a maximal possible accretion 
rate Mmax ~ 31.5MEdd, above which global solutions cannot be constructed as H tends to 
infinity (lines d and e). Here and in Figs. 2, 3, and 6 the vertical hydrostatic equilibrium is 
assumed. 
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Fig. 2. — Advective factor /adv as a function of r. All the lines are the correspondents of 
those in Fig. 1, respectively. Lines d and e are not global solutions as /adv tends to exceed 
its maximal possible value 1. 




Fig. 3. — Total optical depth r as a function of r. Lines a, 6, c, a\ b\ and c' are the 
correspondents of those in Fig. 1, respectively. 
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Fig. 4. — Quantity H/r as a function of r in the case that vertical motion is included. 
Lines a, 6, c, and d are for M/M-^dA — 1, 8.5, 9, and 30, respectively. There is a new 
-^max ~ 8.5MEdd, abovc which there are no global solutions (hues c and d). 




Fig. 5. — Radial velocity Vj. and vertical velocity on the disk surface Vh in global solutions. 
Lines a and b are the correspondents of those in Fig. 4, respectively. 
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Fig. 6. — Specific total energy E as a function of r. Lines a, b, c, d and e are tlie correspon- 
dents of tliose in Fig.l, respectively. 



